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A  Fortran  computer  program  has  been  written  to  compute  the  radar  cross  section 
of  a  conducting  sphere.  The  program  is  useful  when  a  metal  sphere  is  uf°d  as  a  standard 
target  for  calibrating  a  radar  for  target -cross-section  measurement.  It  has  been  used 
to  machine-plot  with  high  precision  a  curve  of  cross  section  (normalized  to  the  optical 
cross  section)  as  a  function  of  the  radius/wavelength  ratio.  The  computed  data  used  to 
plot  the  curve  are  also  presented  in  tabular  form.  The  mathematics  of  the  problem  are 
briefly  reviewed,  and  a  listing  of  the  computed  program  is  given. 
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ABSTRACT 


A  Fortran  computer  program  has  been  written  to  compute  the  radar 
cross  section  of  a  conducting  sphere.  The  program  is  useful  when  a 
metal  sphere  is  used  as  a  standard  target  for  calibrating  a  radar  for 
target- cross- section  measurement.  It  has  been  used  to  machine-plot 
with  hig'i  precision  a  curve  of  cross  section  (normalized  to  the  optical 
cross  section)  as  a  function  of  the  radius /wave length  ratio.  The  computed 
data  used  to  plot  the  curve  are  also  presented  in  tabular  form.  The 
mathematics  of  the  problem  are  briefly  reviewed,  and  a  listing  of  the 
computer  program  is  given. 


AUTHORIZATION 

NRL  Problem  R02-35 
Project  RF-151-402-4011 


CALCULATION  OF  THE  RADAR  CROSS  SECTION  OF  A  PERFECTLY  CONDUCTING  SPHERE 


INTRODUCTION 

This  calculation  is  of  interest  because  spheres  of  high-conductivity 
material  (metal)  are  used  as  standard  targets  for  calibrating  a  radar  for 
target  cross-section  measurement.  Published  curves  based  on  calculations  are 
unsatisfactory  for  precise  work  because  they  appear  to  be  imprecisely  plotted 
(1)  or  the  scale  size  is  too  small  for  accurate  reading  of  values  (2,3). 
Therefore  the  problem  of  writing  a  computer  program  to  calculate  the  cross 
section  and  to  plot  it  accurately  by  machine  was  undertaken. 

The  computation  requires  summing  a  series  whose  terms  involve  spherical 
Bessel  and  Hankel  functions.  These  terms  are  complex  quantities.  Manual 
calculation  would  be  extremely  laborious  and  computer  calculation  also  poses 
some  difficult  problems.  After  a  successful  program  had  been  written,  it 
was  learned  that  others  had  also  dene  it  (3).  Their  programming  approach  was 
somewhat  different,  and  their  actual  program  is  not  presented  in  their  paper. 
The  program  developed  here  is  presented  in  a  form  that  allows  it  to  be.  used 
for  the  calculation  of  cross  sections  of  specific  sphere  sizes,  by  reading 
in  data  giving  the  radius  of  the  sphere  and  the  radar  frequency.  This 
procedure  would  be  used  if  the  curve  and  table  do  not  provide  a  sufficiently 
accurate  number. 

RESULTS 

The  computer  program  was  used  to  machine  plot  a  curve,  Fig.  1,  of  the 
normalized  radar  cross  section  a/ si ra2  ,  where  a  is  the  actual  cross  section  and 
a  is  the  radius  of  the  sphere,  as  a  function  of  the  ratio  a/X,  where  X  is  the 
wavelength.  The  normalizing  quantity  jta2  is  of  course  the  optical  cross 
section  of  the  sphere.  The  plot  was  made  in  a  large  size  (20  x  25  inches)  on 
the  NRL  Gerber  Model  875  Automatic  Drafting  Machine,  with  an  accuracy  of  the 
order  of  .001  inch.  The  coordinate  grid  was  plotted  by  the  machine  also,  so 
that  there  are  no  registration  errors,  as  might  occur  if  the  plot  were  made 
on  standard  graph  paper. 

The  computations  were  made  for  values  of  a/X  from  .05  to  5.045  in  steps 
of  .005.  These  results  are  tabulated  in  Table  1.  The  plot  was  also  made  from 
these  values.  Even  though  the  a/X  interval  .005  is  quite  small,  «.  direct 
digital  plot  does  not  result  in  a  perfectly  smooth  curve.  Therefore  a 
"smoothing"  subroutine,  described  in  a  previous  NRL  report  (4),  was  used  for 
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Table  ij 

i 

Normalized  Values  of  Radar  Cross  Section  (SIGMA)  of  a  Perfectly  Conducting  Sphere 
for  Specified  Radius/Wavel^ngth  (RAT)O)  Values 


MAT  JO 

StSM* 

0.3000 

0.5412 

0.30S0 

0,6544 

0.3100 

0,7606 

0.3150 

0.9158 

0.3200 

1 ,0557 

0.3250 

1.1963 

0.3300 

1.3337 

0.3350 

1.4641 

0.3400 

1  ,5845 

0.3450 

1.6919 

0.3500 

1,7839 

0.3550 

1,8587 

0.3600 

1.9150 

0.3650 

1.9517 

0.3700 

1 .9684 

0.3750 

1.9653 

0.3600 

1.9427 

0.3650 

1.9015 

0.3900 

1.8431 

0.3950 

1.7691 

0.4000 

1 ,6815 

0.4050 

1.5827 

0.4100 

1.4753 

0.4150 

1,3620 

0.4200 

1.2458 

0.4250 

1.1296 

0.4300 

1.0165 

3.4350 

0.9092 

0.4400 

0.8106 

0.4450 

0.7228 

0.4500 

0.6481 

0.4550 

0,5881 

0.4600 

0,5439 

0.4650 

0.5164 

0.4700 

0,5057 

0.4750 

0,5118 

0.4800 

0,5339 

0.4050 

0.5710 

0.4900 

0,6217 

0.4950 

0,6841 

0.5000 

0,7564 

0.5050 

0,8363 

0.5100 

0.9216 

0.5150 

1.0098 

0.5200 

1,0986 

0.5250 

1,1856 

0.5300 

1 .2686 

0.5350 

1.3454 

0.5400 

1.4142 

0.5450 

1.4732 

RATIO 

SIGMA 

0.0500 

0.0860 

0.0550 

0.1254 

0.0603 

0.1767 

0.0650 

0.2418 

0.0700 

0.3229 

0.0750 

0.4218 

0.0800 

0.5404 

0.0850 

0.6800 

0.3900 

0.8417 

0.0950 

1.0258 

0.1000 

1.2315 

0.1050 

1.4572 

0.1100 

1.6999 

0.1150 

1.9549 

0.1200 

2.2164 

0.1250 

2.4772 

0.1300 

2.7294 

0.1350 

2.9646 

0.1400 

3.1747 

0.1450 

3.3524 

0.1S0U 

3.4921 

0.1550 

3.5898 

0.1600 

3.6435 

0.1650 

3.6533 

0.1700 

3.6208 

0.1750 

3.5490 

0.1800 

3.4418 

0.1850 

3.3038 

0.1900 

3.1397 

0.1950 

2.9542 

0.2000 

2.7522 

0.2050 

2.5380 

0.2100 

2.3160 

0.2150 

2.0900 

0.2200 

1 .8638 

0.2250 

1.6412 

0.2300 

1.4255 

0.2350 

1.2202 

0.2400 

1.0284 

0.2450 

0.8533 

0.2500 

0.6976 

0.2550 

0.5641 

0.2600 

0.4548 

0.2650 

0.3718 

0.2700 

0.3161 

0,2750 

0.2884 

0.2800 

0.2886 

0.2850 

0.3159 

0.2900 

0,3688 

0.2950 

— 

0.4449 

RATIO 

SIGMA 

0.8000 

1.1222 

0.0050 

1.0659 

0.8100 

1.0096 

O.MlSo 

0.9550 

0.8200 

0.9034 

0.8250 

0.8561 

0.8300 

0.8143 

0.0350 

0.  .'790 

0.8400 

0.7510 

0.8450 

0.7310 

0.8500 

0.7194 

0.8550 

0.7163 

0.8600 

0.7218 

0.8650 

0.7355 

0.8700 

0.7569 

0.8750 

0.7855 

0.0800 

0.8203 

0.8850 

0.8603 

0.0900 

0.9043 

0.8950 

0.9512 

O.9000 

0.9996 

0.9050 

1.0482 

0.9100 

1.0957 

0.9150 

1.1407 

0.9200 

1.1822 

0.9250 

1.2191 

0.9300 

1.250' 

0.9350 

1.2753 

0,9400 

1.2932 

0.9450 

1.3039 

0.9500 

1.3069 

0.9550 

1.3025 

0.9600 

1.2908 

0.9650 

1.2721 

0.9700 

1.2*7? 

0.9750 

1.2166 

0.9800 

1.1813 

0.9850 

1.1*23 

0.9900 

1.1007 

0.9950 

1.0575 

l.nooo 

1.01*0 

l.no5n 

0.9712 

l.nioo 

0.930* 

1.0150 

0.8925 

1.0200 

0.8586 

1.0250 

0.8294 

1.0300 

0.8057 

1.0350 

0.7881 

1.0400 

0.7770 

1 • 04b0 

0.7725 

i  RATIO 

SIGMA 

^  0.5500 

1.5211 

•  0.5550 

1.5567 

'  0.5600 

1.5793 

:  0.5650 

1.5884 

0.5?00 

1.5840 

■  0.5750 

1.5663 

.  0.5800 

1.5360 

0.5850 

1.49*1 

0.5900 

1.4417 

0.5950 

1.3805 

0.6000 

1.3122 

0.6050 

1.2388 

0.6100 

1.1622 

0.61CJ 

1.0645 

0.&200 

1.0080 

0.6250 

0.934* 

0.6300 

0.8658 

0.6350 

0.0037 

0.6*00 

0.7*98 

0.6*50 

0.7053 

0.6500 

0.6711 

0*6550 

0.6480 

0.6600 

0.6362 

0.6650 

0.6359 

0.6700 

n.6468 

0.6>5o 

0.6685 

0.6800 

0.7001 

0.6850 

0.7407 

0.6900 

0.7009 

0.6950 

0.0435 

0.7000 

0,9028 

0.7050 

0,9653 

0.7100 

1.0290 

0,7150 

1.0O25 

0.7200 

1.1539 

0.7250 

1.2116 

0.7300 

1.2641 

0.7350 

1.3101 

0.7*00 

1.3*83 

0.7*50 

1.3780 

0.7500 

1.3983 

0.7550 

1.4089 

0.7600 

1.4095 

0«765o 

1.400* 

0.7700 

1.3819 

0.7750 

1.3545 

0.7800 

1.3193 

0.78S0 

1.2771 

0.7900 

1.2293 

0.7950 

1 .1772 
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Table  1  (Continued) 


Normalized  Values  of  Radar  Cross  Section  (SIGMA)  of  a  Perfectly  Conducing  Sphere 
for  Specified  Radius/Wavelength  (RATIO)  Values 


RATIO 

SIGMA 

1.0500 

o.rue 

1.0550 

0 

1.0600 

o  ;9«f 

1.0650 

<}.•■.  *5 

1.0TC? 

0.8435 

1.0-50 

0.8758 

1.0000 

0.9096 

1.0050 

0.9*60 

0.98*0 

1.0950 

1.022* 

1.0603 

1.1050 

1.0966 

1.1130 

1.1305 

1.1150 

i.i6:o 

'..12C0 

1.1872 

1 .5  25k 

1.2008 

1  .130.' 

1.2250 

i,l35P 

1  1  2355 

,,1600 

1  • .  ’*C0 

. . ;  *s>' 

i.:  lb 

1.15 

»  .'313 

1.1S'« 

.2183 

1.16 

1.2001 

1  •  lbs ^ 

1.1773 

1.1700 

1.1503 

1.1750 

1.1201 

1.1800 

1.3075 

'.1050 

1.0333 

/ .1900 

1.01H* 

1.1950 

0.9039 

1.200-’ 

0.9506 

1.20b* 

0.9193 

1.2100 

0.8910 

1.2150 

0.8663 

1.2200 

0 .0*58 

1.2250 

0,0301 

1.2300 

0.819* 

1.2350 

0,01*1 

1 .2600 

0.01*1 

1 .2*50 

0.8196 

1.2500 

0.0301 

1.2550 

0.0*55 

1.2600 

0.8652 

1.2650 

0.8807 

1.2700 

0.9152 

1.2750 

0.9**1 

1.2030 

0.97*5 

1.2850 

1.0056 

1.2900 

1.0367 

1.2550 

1.0667 

RATIO 

SIGMA 

1.3000 

1.0951 

1.3050 

1.1209 

1.3100 

1.1*36 

1.3150 

1.1626 

1.3200 

'  .1  77* 

1.3250 

1.1876 

1.3300 

1.1931 

1.3350 

1.1936 

1.3*00 

1.1893 

1.3*50 

1.1803 

1.3500 

1.1669 

1.3550 

1.1*95 

1.360* 

1.1286 

1.3650 

1.10*7 

1.3700 

1.0786 

1.3750 

1.0510 

1.31" ) 

1.0225 

1.3850 

(1.99*0 

1.3900 

(1,9663 

1  3950 

6.9*00 

l.*00O 

(1,9158 

1  .*0bv 

0.89*4 

1.*10  0 

0.8763 

1.*150 

O.8620 

1  .*200 

0.8517 

l.*250 

.0*58 

1  .*300 

0.84*3 

l.*3SQ 

0,8*73 

l.**00 

0.05*5 

1 ,**50 

O.06S9 

1  .*500 

0.0809 

1.4550 

0,8993 

l.*600 

0.920* 

1  .*650 

0.9*37 

1  ,*700 

(1,9685 

l.*750 

0.99*2 

l.*0OO 

1.0200 

1.4850 

1,0*5* 

l.*900 

1.0695 

1.4950 

1.0918 

1.5000 

1.1117 

1.5050 

1.1286 

1.5100 

1.1*23 

1.5150 

1.1522 

1.5200 

1.1582 

1.5250 

1,1602 

1.5300 

1.1681 

1.5350 

1.1520 

1.5*00 

1,1*22 

1.5*50 

1.1289 

RATIO 

Sl«“* 

1.5500 

1.112* 

1.5550 

1.093* 

1.5600 

1.0723 

1.5650 

1.0*96 

1.5700 

1.0260 

1.5750 

1.0022 

1.5800 

0.9T87 

1.5850 

0.9562 

1.5900 

0.9353 

1.5950 

0.9165 

1.6000 

0.9003 

1.6050 

0.8871 

1.6100 

0.8772 

1.6150 

0.8709 

1.6200 

0.8682 

1.6250 

0.8693 

1.6300 

0.87*1 

1.6350 

0.8823 

1.6*00 

0.0938 

1.6*50 

0.9082 

1.6500 

0.9252 

1.655(1 

0.9*42 

1.6600 

1.6650 

0.9861 

1.6700 

1.0079 

1.6750 

1.0295 

1.6800 

1.0503 

1.6850 

1.0697 

1.6900 

1.0*73 

1.6950 

1.1627 

1.7000 

1.1153 

1.7050 

1.12*9 

1.77,0 

1.1313 

1.7150 

1.13*3 

1.7200 

1.1338 

1.7250 

1.1300 

1.7300 

1.1229 

1.7350 

1.1128 

1.7*00 

1  .0999 

1.7*50 

1.08*6 

1.7500 

1.067* 

1.7550 

1  .0*86 

1.7600 

1.0289 

1.7650 

1.0007 

1.7700 

0.9887 

1.7750 

0.9692 

1.7800 

0.9509 

1.785c 

0.93*2 

1.7900 

0.9195 

1 . 7950 

1 _ 

0.9072 

RATIO 

SIGMA 

1.0000 

0.8977 

1.0050 

0.8911 

1.0100 

0.8876 

1.0150 

0.8873 

1.0200 

0.8901 

1  •  *250 

0.8959 

l.«300 

0.90*7 

1.8350 

0.9160 

1.8*00 

0.9297 

1.8*50 

0.9*52 

1.0500 

0.9622 

1.0550 

0.9802 

1.0600 

0.9980 

1.0650 

1.0173 

1.0700 

1.035* 

1.0750 

1.0525 

1.0800 

1.0683 

l.«850 

1.0822 

1.8900 

1.09*0 

1.0950 

1.1032 

1.9000 

1.1099 

1.9050 

1.1136 

1.9100 

1.11*5 

1.9150 

1.112* 

1.9200 

1.1075 

1.9250 

1.0998 

1,9300 

1.0897 

1.9350 

1.077* 

1.9*00 

1.0633 

1.9*50 

1.0*77 

1.9500 

1.0311 

1.9550 

1.01*0 

1,9600 

0.9966 

1.9650 

0.9797 

1.4700 

0.9635 

1.9750 

0.9*85 

1.9800 

0.9352 

1,9860 

0.9237 

1.9900 

0.91*5 

1  .9950 

0.9078 

2.0000 

0.9037 

2.(1050 

0.9022 

?>0100 

0.9035 

2.0150 

0.907S 

2,(1200 

0.9140 

2.0250 

0.9229 

2. 0300 

0.9338 

2,0350 

0.9*66 

2.0*00 

0.9607 

2,O*S0 

0.9760 

Table  1  (Continued) 
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Table  1  (Continued) 

Normalized  Values  of  Radar  Cross  Section  (SIGMA)  of  a  Perfectly  Conducting  Sphere 
for  Specified  Radius/Wavelength  (RATIO)  Values 


RATIO 

SIGMA 

3.S500 

0.9620 

3.5550 

0.9636 

3.5600 

0.9662 

3.5650 

0.9697 

3.5/00 

0.97*0 

3.5750 

0.9789 

3.5600 

0.90*5 

3.5950 

0.990* 

3.5900 

0.9966 

3.5950 

1  .0029 

3.6000 

1.0090 

3.6050 

1.0150 

3.6100 

1 .0205 

3.6150 

1.0255  , 

3.6200 

1.6298 

3.6250 

1.0333 

3.6300 

1.0359 

3.6350 

1.0376 

3.6*00 

1  .036* 

3.6*50 

1.0181 

3.6500 

1.0368 

3.6550 

1.03*6 

3.6600 

1.0315 

3.6650 

1.0276 

3.6/00 

1.0230 

3.6750 

1.0179 

3.6600 

1.0123 

3.6650 

1.0065 

3.6900 

1.0005 

3.6950 

0.99*6 

3. 7000 

0,9*89 

3.7050 

0.9*35 

3.7100 

0.9786 

3.7150 

0.97*3 

3.7200 

0.9707 

3.7250 

O.9680 

3.7300 

0.9661 

3.7350 

0.9651 

3.7*00 

0.9651 

3.7*50 

0.9661 

3.7500 

0.9679 

3.7550 

0.9706 

3.7600 

0.97*t 

3.7650 

0,9783 

3.7700 

0.9*31 

3.7750 

0,9883 

3.7600 

0,9938 

3.7650 

6.9*95 

3. 7V00 

1 .0051 

3.7950 

1 .0106 

RATIO 

SIGH  A 

3.8000 

1.0159 

3.«050 

1.0207 

3.8100 

1.02*9 

3.8150 

1.0285 

3.8200 

1.031* 

3.8250 

1.033* 

3.8300 

1.03*6 

3.«350 

1.03*8 

3.8*00 

1.03*1 

3.8*50 

1.0326 

3.8500 

1.0302 

3.8550 

1.0271 

3.8600 

1.0233 

3.8650 

1.0189 

3.8700 

1.01*0 

3.8750 

1.0088 

3.8800 

1.003* 

3.8850 

0.9980 

3.8900 

0.9927 

3.8950 

0.9876 

3.9000 

0.9829 

3.9050 

0.9787 

3.9100 

0.9750 

3.9150 

0.9721 

3.9200 

0.9700 

3.9250 

0,968* 

3.93U0 

0.9681 

3.9350 

0.5685 

3.9*00 

0.9698 

3.9*50 

0.9718 

3.9500 

0.97*6 

3.9550 

0.9781 

3.9600 

C.9B21 

3.9650 

0.9866 

3.9700 

C.9915 

3.9750 

0.9966 

3.9800 

1.0018 

3.9850 

1.0069 

3.9900 

1.0118 

3.9950 

1.016* 

*,0000 

1.0206 

*.O0S0 

1.02*2 

*.0100 

1.0272 

*.0150 

1.0295 

*.0200 

1.0310 

*.0250 

1.0316 

*.0300 

1.0315 

*.0350 

1.0306 

*.o*oo 

1.0288 

*,0«50 

1.0263 

HA  1 10 

. .  '  1 

STGHA 

3.3000 

1.0075 

3.3050 

1 .0002 

3.3100 

".9930 

3.3150 

n.9860 

3.3200 

ft,979* 

3.3250 

6.9735 

3.3:10 

6.9683 

3.3350 

6,96*0 

3.3*00 

(1.9607 

3.3*50 

6,9585 

3.3500 

6.95/5 

3.3550 

6.9576  i 

3.360b 

6.9588 

3.3650 

6.9612 

3.3700 

6.96*6  I 

3.3750 

6.9689 

3.3800 

6.97*1 

3.3650 

6.9800 

3.3900 

6.986* 

3.3950 

6.9931 

3. *000 

1 .0000 

3. *050 

1.0069 

3. *100 

1.0136 

3. *150 

1.0200 

3. *200 

1.0258 

3. *250 

1.0300 

3. *300 

1.0352 

3. *350 

1  .0385 

3. **00 

1.0*09 

3. **50 

1.0*22 

3. *500 

1  .0*2* 

3. *550 

1.0*15 

3. *600 

1.0395 

3. *650 

1 .0365 

3. *700 

1 .0326 

3. *750 

i  .0279 

3. *600 

1  ,0225 

3. *650 

1.0165 

3. *900 

1,0102 

3. *950 

1 .0036 

3.5000 

6.9971 

3,5050 

6.9906 

3.5100 

6.98*5 

3.5150 

6,9788 

3.5200 

6.9737 

3,5250 

6,969* 

3.5300 

6,9660 

3.5350 

6.963* 

3.5*00 

6.9619 

3.5*50 

6.9615 

RATIO 

SIGMA 

3.0500 

1.0*79 

3.0550 

1.0508 

3.060O 

1.0523 

3.0650 

1.052* 

3. 070W 

1.0513 

3.0750 

1.0*87 

3.0800 

1.0*50 

3.0850 

1.0*01 

3.090O 

1.03*2 

3. 095O 

1.0275 

3.1000 

1.0201 

3.1050 

1.0123 

3.1100 

1.00*2 

3.1150 

0.9961 

3.1200 

0.9881 

3.1250 

0.9806 

3.1300 

0.9737 

3.1350 

0.9675 

3.1*00 

0.9622 

3. l*So 

0.9580 

3.150O 

0.9550 

3.1550 

0.9532 

3. 1600 

0.9527 

3.1650 

0.9535 

3.1700 

0.9556 

3.1750 

0.9588 

3.1800 

0.9632 

3.1850 

0.9685 

3.1900 

0.97*7 

3.1950 

0.9815 

3.2000 

0.9888 

3.2050 

0.996* 

3.2100 

1.00*1 

3.2150 

1.0117 

3. 220o 

1.0189 

3.2250 

1.025/ 

3.2300 

1.0317 

3.2350 

1.0370 

3.2*00 

1.0*12 

3,2*So 

1.0**3 

3.250O 

1.0*63 

3.25SO 

1.0*71 

3.2600 

1.0*67 

3.2650 

1.0*50 

3.2700 

1.0*22 

7.2750 

U038* 

3.2800 

1.0335 

3.28SO 

1.0279 

3.2900 

1.0215 

3.2950 

1.01*7 

Table  1  (Continued) 


Normalized  Values  of  Radar  Cross  Section  (SIGMA)  of  a  Perfectly  Conducting  Sphere 
for  Specified  Radius/Wavelength  (RATIO)  Values 


RATIO 

SIGH  A 

*.•<000 

1.0223 

*.«030 

1.0219 

*.**100 

1.0208 

*.**130 

1.0192 

*.**?oo 

1.0172 

*.8250 

1.01*7 

*.*300 

1.011** 

*.•*330 

1.0086 

*.*»*0o 

1.0053 

♦  .*>*30 

1.0018 

*.*>500 

0.9383 

*.•<530 

0.99*9 

♦  .**600 

0.9916 

*.**630 

0.9e86 

*.**700 

0.9859 

*.*>730 

0.9836 

*.*>H00 

0.9817 

*.*>**30 

0.980* 

*.•*900 

0.9796 

*.*>930 

0.9793 

*.9000 

0.9796 

*.9030 

0.980* 

*.9100 

0.9818 

*.9130 

0.9836 

*.9200 

0.9859 

*.9230 

0.9686 

*.9300 

0.9915 

*.9350 

0.99*7 

*.9*00 

0.998c 

♦.9*50 

1.001* 

*.9500 

1.00*7 

♦.9530 

1.0079 

*.9600 

1.0109 

*.9630 

1.0136 

*.9700 

1.0159 

*,9730 

1.0178 

*.9000 

1.0192 

*.9830 

J .0202 

*.9900 

1.0206 

*,9930 

1.0206 

5.0000 

1.C198 

5.0030 

1.0187 

5.O100 

1.0170 

5.0130 

1.01*9 

5.0200 

1.0125 

5.0230 

1.0097 

5.0300 

1.0067 

5.0330 

1.0036 

5.0*00 

1.000* 

5.0*30 

0.9971 

RATIO 

SlOHA 

*.5500 

0.99*5 

*.5550 

0,998* 

*.S600 

1.002* 

*.5650 

1.0063 

*.5  TOO 

1  .0100 

*.5730 

1  .0135 

*.5600 

1 .0166 

*.5830 

1 .0192 

*.5900 

1.021* 

*.5930 

1.0230 

*•6000 

1.0239 

*•6050 

1.02*3 

*•6100 

1 .o?*0 

*.6130 

1.0232 

*•6200 

1.0217 

*.6230 

1.0196 

*.6300 

1.0171 

*.6350 

1.01*1 

*.6*00 

1*0108 

*•6*50 

1 .0072 

*•6300 

1 .0035 

*•6350 

0.9997 

*.6600 

0.9959 

4*66$0 

0.9923 

*.6700 

0.9*89 

*.67so 

0.9*58 

*.6800 

0.9831 

*.685n 

0,9*09 

*.6900 

0.9792 

*.6950 

0.9781 

*.7000 

0.9776 

*.7030 

0.9776 

*.7100 

0.9763 

*.7150 

0.9796 

*•7200 

0.981* 

*. 723C 

0.983' 

*.7300 

0.9*6* 

♦.7330 

0,9*95 

*.7*00 

0,9928 

*.7*S0 

0.996* 

*.7500 

1.0000 

*.7550 

1.0037 

*.7600 

1 .0072 

*.7630 

1 .0105 

*.7700 

1.01 76 

*.77sn 

1 .0163 

*.7600 

1 .0185 

*.7850 

» .0203 

*,79on 

1.0216 

*.7yso 

1.0222 

RATIO 

MSHA 

*.3000 

6.9770 

*.3050 

•>.9751 

*.3100 

e.9739 

*.3150 

0,973* 

*.3200 

0,9736 

*.3250 

0.97*5 

*.3300 

0.9761 

*.3350 

0.978* 

i  *.3*00 

0.9812 

|  *.3*50 

0.98*5 

*.3500 

0.9882 

*.3550 

0.9923 

*.3600 

0.9965 

*.3650 

1.0008 

*.3700 

1.0051 

*.3750 

1.0093  | 

*.3800 

1.0132  | 

*.3*S0 

1.0167 

*.3900 

1.0198 

*•3950 

1  .022* 

*.*000 

1.02** 

*.*050 

•.0257 

*•♦’00 

1 .026* 

*•*150 

1  .026* 

*•*200 

1.0257 

*.*250 

1.02** 

*.*300 

1.022* 

*.♦350 

1.0198 

*.**00 

1.0166 

*.**50 

1,0133 

*.*500 

1  .0095 

*.*550 

1.0055 

*.*600 

1.0C1* 

*.*650 

0.9973 

*.*700 

fl.9932 

*.♦750 

0,989* 

*.*■00 

0,9859 

*,*■50 

0.9828 

*.*900 

0,9802 

*.*950 

0.9781 

*.5000 

fl.9766 

*.5(l50 

0.9758 

*.5100 

n.9756 

*.5150 

n.9761 

*.5200 

n.9772 

*.5250 

0.9789 

*.5300 

n.9812 

*.5350 

fl,98*0 

*•5*00 

0.9872 

*.5*5o 

0,9907 

R»TIO 

sigma 

*.0500 

1.0232 

*.0550 

1.0195 

4,0600 

1.0153 

*.0650 

1.0107 

*.0700 

1.0059 

*.0  750 

1.0010 

*•(1800 

0,9960 

*.0850 

0,9913 

*.090u 

0.9867 

*.095o 

0.9826 

*.1000 

0.9790 

*.1050 

0.9760 

*.1100 

0.9736 

*.1150 

0.9720 

*.l20u 

0.9711 

*.1250 

0.9710 

*.1300 

0.9717 

*.1350 

0.9731 

*.1*00 

0.9753 

*,l*5o 

0.9781 

6.l5oo 

0.9815 

6.1550 

0.985* 

6.1600 

0.9897 

♦.1650 

0.99*2 

*.1700 

0.9989 

*.1750 

1.0036 

*.1«00 

1.008? 

*.1850 

1.0126 

*.1900 

1.0167 

*.1950 

1.0203 

*.2000 

1.0233 

*.2050 

i .0258 

♦.2100 

1.0276 

*.2150 

1.0286 

*.2200 

1.0289 

*.2250 

1.0285 

4.2300 

1.0273 

*.2350 

1.025* 

*.2*00 

1.0229 

*.2*50 

1  .C 198 

*.2500 

1.0162 

4.2550 

1.0122 

4.2600 

1.0079 

*.2650 

1.003* 

*.2700 

0.9989 

♦.2750 

0.99*5 

*.2«00 

0.9902 

4.2850 

0.9862 

*.2900 

0.9826 

*.?95o 

0.9795 

the  plotting;  this  subroutine  interpolates  additional  points  where  they  are 
needed,  using  a  coM.c  interpolating  polynomial. 

The  program  was  then  modified  so  that  it  calculates  the  radar  cross 
section,  in  both  square  meters  and  square  feet,  of  a  specific  sphere  at  a 
specific  frequency.  Printouts  of  sample  calculations  using  this  program  are 
shown  as  Figs.  2,  3,  and  4.  A  listing  of  the  program  is  given  later  in  this 
report. 

MATHEMATICAL  FORMULAS 

Kerr1  gives  the  following  equation  (p.  451)  for  the  radar  cross  section 


2 

(1) 


a  of  a  sphere  of  radius  a: 


a  _  1 

— S'  "  T 
jta  p 


OO 

£  M)"  (2n  +  1)  0»nS  "  b„8  ) 


n=i 


where  p  =  2na/X,  X  is  the  wavelength,  and  the  a^  and  bns  are  terms  of  a 
"multipole  expansion"  —  that  is,  these  terms  are  proportional  to  the  ampli¬ 
tudes  of  magnetic  and  electric  multipoles  induced  in  the  sphere  by  the  inci¬ 
dent  wave.  When  the  sphere  is  perfectly  conducting  (Kerr,  p.  452): 


s 

% 


in  (P) 
hn(2)  (P) 


-  C  P  Jn  (P)  ]’ 

[  P  hn<s>  (p)  ]• 


(2) 

(3) 


where  the  primes  denote  differentiation  with  respect  to  the  argument.  The 
functions  jn  and  h^2)  are, respectively ,  the  spherical  Bessel  function  of 
the  first  kind,  and  the  spherical  Hankel  function  of  the  second  kind. 


A  subroutine  to  evaluate  jn,  jn',  and  also  the  spherical  Bessel 
function  of  the  second  kind  and  its  derivative,  yn  and  yn' ,  was  obtained 
from  the  NRL  CDC-3800  program  library.*  No  subroutine  for  evaluating  the 


*Library  Catalog  Identification  C3-UCSD-BFFGH.  Subroutine  BFFGH  was  written 
by  Frank  Hagin  of  Texas  Instruments,  Dallas  Texas. 
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Li ZaaL*.  ■  m  waii: 


2 


CAl.CU.AT J0N  0F  RADAR  CROSS  SECTIEN  OF  SPHERE 
REFFRENCE  -  KERR i  *PRflP AG AT ! £N  OF  SH0RT  RAD10  WAVES* 
RAD  LA8  SERIES  V0L ,  13,  P,  451 ,  EQ,  29, 


INPUT  DATA  AND  DERIVED  QUANTITIES  -- 


RADIUS,  METERS  1,5?4 

RADIUS,  FEET  . . 5,000 

RADAR  FREQUENCY,  MEGAHERTZ  . .  138,60 

WAVELENGTH,  METERS  , . ?,16 

RADIUS/WAVELENGTH  RATIP  ,  ,  ,  0,7046 

CJKCUMFERENCfc/WAVELEA'GTH  RATIO  .  4,4270 


eutput  DATA  — 


NUMBER  0F  SERIES  TER^S  ADCFC  .  12 

RADAR  CROSS  SECT  1 0N,  NORMAL I?ED  T0  OPTICAL  0,9599 

RADAF  CROSS  SECTION,  SQUARE  MfcTERS  .  7,0038 

RADAR  CROSS  SECTION,  SOlARF  FEFT  .  75,3884 


Fig.  2  -  Computer  printout  from  Program  SPHERE  for  radius 
5  feet  and  frequency  138.6  MHz. 


CALCULATJfN  £F  "ADAH  CROSS  S£r 7 1 tN  0F  SPHERF 
REFERENCE  -  HERR,  *PRWPAGaTIGN  OF  SHWRT  RADIO  WAVES* 
RAO  LAW  SfRlfcS  V«L,  13,  P,  451,  fcfJ,  29, 


INPUT  DATA  A!,0  PcRlVtn  QUANTITIES  -- 


RADIUS,  CENTIMETERS  . 5)0,00 

RADIUS,  FEET  . . 1.640 

radaf  frequency,  megahertz  .  loon.ou 

WAVELENGTH.  CENTIMETERS  . . 29,98 

RADIUS/WAVF LENGTH  RATin  . .  1,6678 

ClRCUMFtRE'lCE/WAVfcLENGTh  RATIO  .............  10,4792 


output  data  -- 

numheh  sf  series  terms  ahcfl  .  ?o 

RADAR  CROSS  SECTION,  NORMALIZED  Ta  OPTICAL  0,9904 

RADAR  CROSS  SECTION,  SQUARE  METERS  .  0,7841 

SCU4MF  MfcTERS,  EXPONENT  FORMAT  . ,,,,7.8411-001 

RADAR  CROSS  SECTION,  SQUARE  FEFT  . .  8,4401 


Fig.  3  -  Computer  printout  from  Program  SPHERE  for  radius 
0.5  meter  and  frequency  1000  MHz. 


calculation  er  radar  cross  sectisn  of  sphere 

REFERENCE  *  KERR,  ^PROPAGATION  OF  SHORT  RADIO  WAVES* 
RAD  l,Ag  SERIES  VOL,  1 3,  P,  451,  EQ,  29, 


INPUT  DATA  AND  DERIVED  QUANTITIES  -- 


RADIUS,  CENTIMETERS  . . 

RADIUS,  INCHES  . 

RADAR  FREQUENCY,  MEGAHERTZ 
WAVELENGTH,  CENTIMETERS 
RADIUS/WAVELENGTH  RATIO 

circumference/wavelength  RAT Iq 


6,00 
2,  36 

1000,00 

2R.98 

0,2001 

1,2575 


output  DATA 


number  of  series  TERMS  AT. DEC  . . .  6 

RADAR  CROSS  SECTION,  NORMALIZED  Ta  OPTICAL  2,7464 

RADAR  CROSS  SECTION,  SQUARE  METERS  .  0,0311 

SQUARE  METERS,  EXPONENT  FORMAT  3, 1061-002 

RADAR  CROSS  SECTION,  SQUARF  FfcET  0,3343 

SQUARE  FEET,  EXPONENT  FORMAT  . 3,3434-001 


Fig.  4  -  Computer  printout  from  Program  SPHERE  for  radius 
0.06  meter  and  frequency  1000  MHz. 
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spherical  Hankel  function  was  available.  However,  the  problem  was  solved 
by  the  use  of  an  expression  relating  the  spherical  Hankel  function  to  the 
spherical  Bessel  functions.  The  necessary  relation  is  given  by  Messiah  (5) 


hk(2)  *  ik  '  ^k 

hk(2)’  *  V  -  ^k' 

Since  Subroutine  BFFGK  gives  jfc,  ,  y^>  and  y^' ,  it  is  a  simple  matter  to 
obtain  Uk^3'  and  hk^2'  •  (To  interpret  the  relations  given  by  Messiah,  it 
is  necessary^fcc  note  that  his  function  nk  is  the  negative  of  yk»  and  his 
functions  hk',+^  and  '  are  equal  to  ihk^1^  and  -ihk'2^,  respectively.) 

For  very  small  values  of  p,  the  approximation* 


-A-  =  9p4  =  1.4027  (a/X)4  x  104  (4' 

ita 

can  be  used;  this  is  the  well  known  Rayleigh  scattering  law  (for  back 
scattering) . 

For  very  large  values  of  p  (or  a/X) ,  the  asymptotic  result  is 


1 


As  will  be  discussed  in  some  detail,  this  approximation  is  valid  to  4  signi¬ 
ficant  figures  for  a/X  >  44  (  and  possibly  for  somewhat  smaller  values). 

FORTRAN  PROGRAM 


The  program  written  to  compute  the  cross  section  is  simply  a  Fortran 
algorithm  to  evaluate  Eq.  (1).  The  summation  of  multipole  terms  (which  are 
complex  numbers)  is  carried  to  the  point  to  which  the  fractional  change  in  a 


*Reference  1,  page  452 


A1 '  w<  3J  **v  yci"  ■ 


due  to  the  last  terra  added  is  less  than  10 ~6 .  No  exact  analysis  of  the 
accuracy  obtained  by  this  procedure  was  made;  however,  a  check  on  accuracy 
is  afforded  by  the  knowledge  that  as  a/X  becomes  large,  the  envelope  of  the 
oscillations  of  a  asymptotically  approaches  jtaa.  The  numerical  results 
obtained  for  large  a/X  agree  with  this  prediction.  To  test  this,  the  program 
was  run  with  the  value  of  a/X  doubled  for  each  successive  calculation.  The 
results  are  shown  in  Table  2.  The  initial  value  of  a/X  was  .0850  and  this 
was  doubled  until  the  value  21.76  was  reached,  at  which  time  the  run  terminated 
An  additional  run  was  then  made  for  the  values  a/X  «  44  and  a/X  -  88.  The 
second  column  of  the  table  gives  the  number  of  series  terms  summed  before  the 
above-stated  criterion  was  met  (fractional  change  in  a  <  10"6). 


Table  2 


Results  for  Large  Values  of  a/X 


a/X 

No.  of  Terms  Summed  , 

a/jta3 

.0850 

5 

0.6800 

.1700 

6 

3.6208 

.3400 

7 

1.5845 

.6800 

11 

0.7001 

1.3600 

17 

1.1286 

2.7200 

28 

1.0156 

5.4400 

47 

0.9929 

10.8800 

83 

0.9970 

21.7600 

158 

0.9997 

44.0000 

303 

1.0000 

88.0000 

586 

1.0000 
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It  is  seen  that  for  a/X  5  44,  the  values  of  cr/jta2  differ  from  1  by 
less  than  5  x  10*6 .  This  suggests  that  the  calculation  is  accurate  to  at 
least  4  significant  figures  for  values  of  a/X  up  to  at  least  88.  To  make 
sure  that  the  values  1.0000  were  not  obtained  "by  accident"  at  a/X  *  44, 
the  program  was  run  with  values  a/X  *  44.04  and  44.045.  These  numbers  were 
chosen  to  insure  that  if  a/X  just  happened  to  result  in  a/ ita2  =  1.0000  when 
the  value  is  still  actually  oscillating  appreciably;  the  oscillation  would 
be  revealed.  The  results,  printed  out  to  additional  decimal  places,  were 
a/ira2  *=  0.999982  and  0.999981,  respectively.  It  thus  seems  reasonable  to 
conclude  that  the  calculation  does  have  at  least  4-significant-figure 
accuracy. 

Since  for  large  values  of  a/X  the  number  of  terms  summed  almost  doubles 
with  each  doubling  of  a/X,  a  point  would  eventually  be  reached  at  which  the 
accuracy  would  be  reduced  by  computer  round-off  error.  However,  in  the 
CDC-3800  this  event  occurs  at  a  value  of  a/X  well  above  the  point  at  which 
the  approximation  cr/jta2  «  1  can  be  assumed  to  4  significant  figures.  It  is 
therefore  recommended  that  for  values  of  a/X  >  44,  the  optical  approximation 
cr/jta2  =  1.000  be  used  instead  of  making  an  actual  computer  calculation. 

The  time  required  for  the  calculation  with  the  CDC-3800  is  approximately 
58  milliseconds  per  series  term,  not  including  compilation  time  and  input- 
output  operations.  Therefore  calculation  of  single  values  of  a  for  the 
usual  range  of  values  of  a/X  requires  only  a  few  seconds  at  most. 

A  listing  of  the  Fortran  program  that  has  output  of  the  type  shown  in 
Figs.  2-4  follows.  The  input  to  this  program  is  a  data  card  or  cards  (one 
for  each  combination  of  sphere  size  and  frequency  to  be  calculated) .  Any 
number  of  data  cards  can  be  stacked  at  the  appropriate  place  in  the  deck, 
and  followed  by  an  end-of-file  card.  The  program  will  then  read  each  card, 
do  the  appropriate  computations,  and  print  out  the  results.  The  data  card 


format  is: 

Card  Columns 

Quantity 

Format  Specification 

1-10 

Radius  of  sphere 

F10. 

11-20 

Radar  frequency 

F10. 

30 

Radius  units  designator 

11 

14 


The  radius  of  the  sphere  can  be  given  in  either  meters  or  feet.  If  it 
is  given  in  meters,  Column  30  is  left  blank.  If  it  is  given  in  feet,  a  1  is 
punched  in  Column  30.  The  radar  frequency  is  given  in  megahertz.  Both  the 
radius  and  the  frequency  can  be  punched  as  floating-point  numbers  anywhere 
within  the  specified  fields;  a  decimal  point  must  be  included  unless  the 
numbers  are  integers  ar.d  are  right-adjusted  in  the  field. 

As  mentioned  in  the  introduction,  the  computer  program  that  follows  was 
written  before  the  paper  of  Adler  and  Johnson  (3)  was  called  to  the  author’s 
attention.  They  describe  a  method  of  computing  the  spherical  Bessel  and 
Hankel  functions  based  on  a  recursion  relation.  Use  of  this  method  instead 
of  Subroutine  BFFGH  would  have  been  somewhat  more  direct  and  possibly  would 
reduce  computing  time;  however,  since  the  program  as  it  stands  works  well  and 
runs  without  excessive  use  of  computer  time,  it  was  decided  not  to  change  it. 
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PROGRAM  SPHERE 

COMPUTES  RADAR  CROSS  SECTION  OF  A  PERFECTLY  CONDUCTING  SPHERE, 
MATHEMATICS  GIVEN  BY  KERR,  *PR0PAGATI§N  OF  SHORT  RADIO  WAVFS* 

RAO  LAO  SERIES  VOU,  13  (MCGRAW-HILL,  1950),  P.  451,  EO,  29, 

PROGRAM  WHITTEN  BY  l,,  V.  BLAKE,  NRL  CODE  5370,  THIS  VERSION  3/1/72 

COMPUTES  AS  MANY  CASES  (SPHERE  SUES,  FREQUENCIES)  AS  THERE  ARF 
DATA  CARDS,  RADIUS  OF  SPHERE,  COLS,  1-10 ,  FREQUENCY  MEGAHERTZ, 
COLS.  11*20 ,  IF  RADIUS  IS  GIVEN  IN  FEET,  PUNCH  A  1  IN  COL,  30, 

If  RADIUS  IS  GIVEN  JN  MfcTfcRS,  NO  PUNCH  (OR  ZERO  PUNCH)  IN  COL.  30, 

TYPE  COMPLEX  HN2,UHN2#ANS,0NS, TERM, SUM 
DATA  (PI*3,  1415926536),  (MiO) 

500  READ  30,  RADIUS,  f  RQNCY,  JUNJTS 
IF  (EOF ,60 )  80,89 
08  STOP 

09  If  UUNITS,EG,1>  31,32 

31  RADFT  «RAD jUS 
RADIUS*RAUIUS«.3048 
GO  TO  33 

32  RADFT«RADIUS/,3Q40 

33  WVLGTH«29V,7925/FHQNCY 
RATIO*  RADIUS/WVLUTH 
RHO  *  2 ,  *P  I *RAT 1 0 

N  s  0 
QSQ  s  1, 

IKO  s  1 
Y  «  1, 

N  *  0 

SUMLAST  *  0, 

ISIGN  *  1 

SUM  s  CMPLX ( 0 , , 0 , ) 

IN*  N*1 

ISIGN  «  IS!GN*(«1) 

CALL  BFFGH(RH0,N,BJ,DdJ,DURj,OY,DBY,DDBY,QSQ, IQK, Y) 

HN2  *  CMPLX ( 8 J, *0 Y  ) 

DMN2  ■  CMPLX(DBU,*nBY) 

ANS  =  -BJ/HN2 

BNS  *  *(RhO*DBJ*8U)/( HHy*UHN2*HN2 ) 

TERM  *  ISIGN* (2*N*1)*(*NS-BNS) 

SUM  s  SUM*TERM 
ASSUME  s  (CABS(SUM) )••? 

DlFFRsARS  ( (A8SUM<!*SUMLAST>/AbSUM2) 

IF  (DIFFR  ,UT,  l,t-6)  GO  TO  9 
SUMLAST  *  A8SUM2 
GO  TO  1 

9  SJGNML  «  ABSUM2/(RHO*RHO) 

SIGMA*aIGNML*Pl*RADIUS*RAUIUS 
SIGFT*SIGMA*10, 76391 
If  (M | EQ . 0 )  80,81 
eo  m*i 

GO  TO  82 
01  PRINT  998 
0?  PRINT  992 
PRINT  99 
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PRINT  990 
PRINT  991 
PRINT  999 
PRINT  100 

If  (RADIUS  ,GE.  1,)  4U ,  41 

40  PRINT  101 «  R*U I U? 

43  PRINT  102,  RADPT 
GO  T0  42 

41  RADCM«RADIUS«100. 

PRINT  1011,  RADCM 

If  (RADFT  ,GE .  1.)  43,44 

44  RADJN»RADFT*12. 

PRINT  1012,  RAOJN 

42  PRINT  103 | f RQNCy 

If  (WVLGTH.GE.l , )  5Q,i>l 

50  PRINT  104,  wvugth 
G0  TO  52 

51  UVLGHH"WVIGTH*100. 

PRINT  105, WVUGHH 

52  PRINT  106,  RAT  1 0 
PRINT  107,  RM0 
PRINT  108 
PRINT  1080,  N 
PRINT  109,  SIGNMU 
PRINT  110,  SIGMA 

If  (SJGMA  ,11,  1,  , 0R ,  SIGMA  ,GT ,  99999,9999)  PRINT  11?,  SIGMA 
PRINT  111,  SIGFT 

If  CSJQFT  ,UT,  1,  ,0R,  SIGFT  ,QTf  99999,9999)  PRINT  113,  SIGFT 
G0  T0  500 

40  F0RMAT<2F1O,2,9X,I1) 

998  F0RMAT(1M1) 

992  FORMAT (//////) 

99  FORMAT (19X, ‘CALCULATION  0)  RADAR  CR0SS  SECTION  0F  SPHERE*  > 

990  F0RMAT(15X,*RfcFfcRfcNCF.  -  KERR,  PPR0PAGAT  J0N  0F  SHORT  RADI®  WAVES/* ) 

991  FORMAT (15X,*RAD  J.AB  SERIES  VOL ,  13,  P,  451,  FQ,  29,  •//) 

999  F0RMATU5X,* . - . •///> 

100  F0RMAT(15X, »INPUT  DATA  AND  DERIVED  QUANTITIES  --*/> 

101  F0RMAT(2OX,*RAD?US,  METERS  . . . . . .  ,,,,*,  Fin,  3  ) 

1011  FORMAT (20X , *RADI US ,  CENTIMETERS  . . ,..*,F10.?  ) 

1012  FORMAT  (20X ,  *RADJ  US  ,  INCHES  . . .  ) 

10?  FORMAT  (20X,*RA0JUS,  FbFT  . . ,,.*,F10,3  ) 

103  FORMAT(20X,*RADAR  frfuuency,  megahertz  . . *,nn,?  ) 

104  FORMAT  (20X,  *  WAVE  LENGTH,  METERS  . . ,.*,F10,?  ) 

105  FORMAT (20X,*WAVELfcNQTH,  CENTIMETERS  ,»,F10,2  ) 

106  FORMAT (20X , *RADIUS/4AVELENGTH  RATIO  . . . Flf) , 4  > 

107  FORMAT<20X,*CIRCUMF6RENGF/WAVELENGTH  RATIO  . . . . . . *F10 , 4/// ) 

108  FORMAT <15X,*0UTPUT  DATA  --•/> 

1080  F0RMAT(2CX,*NUMBEH  OF  SERIES  TERMS  ADDED  . . . 5X , 1 5  ) 

109  F0RMAT(20X,*RADAR  CRSSS  SECTION,  NORMALIZED  TO  OPTICAL  *,F10,4  ) 

110  FORMAT (20X,*RADAR  CRASS  SECTION,  SQUARE  METERS  . t,.*,F10,4  > 

112  FORMAT (23X,*SQUARE  MEIERS,  EXPONENT  FORMAT  . ,.,*,£10,4  ) 

111  FORMAT (2 OX, ‘RADAR  CROSS  SECTION,  SQUARE  FEET  . ,*,F10,4  > 

113  FORMAT (23X,*SUl>ARb  FEET,  EXPONENT  FORMAT  . ,*,E10,4  ) 

E  NO 


17 


03/06/72 


SUBROUTINE  BFFGH  (Z,  L,  BA,  BjP,  BJPP,  BB,  BYP,  BYPP,QSQ, JKQ, Y 
C  C3  UC$0  BFFQH 

D1MENSI0N  BJ(b>,  BY(5),  BYL(5),R<3) 

YSAV6  *  Y 

If  <Y  •  1,)  102,  102,  101 

101  Y  *  1, 

102  FL  *  L 
Ll«l*i 
Z8  a  Z 

n  «  o. 

FLl  *  FL  *  1. 

ZSQ  a  Z  •  Z 

FLR  a  SQRTF  (FL  •  FLl) 

SINE  ?  SINF(Z) 

cbse  »  cssf(7) 

IF  ( IKQ  -  1)  104,104,10b 
100  IF  <Z  -  FLR)  800,«00,VftO 
104  IF  (Z«l,)  300,200,200 

109  IF  (QSQ)  110,  110,  104 

110  IF  (Z  •  1,)  500,  400,  400 
200  80(1)  *  SINE  /  ZB 

SIGN  «  1. 

ER  a  -Y 

BY ( 1 )  *  ER  •  C0SE 

210  BY  ( 1 )  a  BYU)  /  Zb 
ER  a  «J(1)  /  ZB 

IF  (SIGN)  211,212,212 

211  ER  >  .ER 

212  8J<2>  ■<8Y<1)/Y  ♦  ER)/  Y 

BY  ( 2 )  *  (dY(l)/Z  -  Y*dJ(l))*Y 
FLX  *  2, 

221  TtMP  f(2,  *FLX  -1,)/Z 
BJ(3)1(TEMP  •  BJ(2i  -  BJ(1>/Y)/Y 
IF  (SIGN)  223,224,224 

223  BJ(3)  «  -BJ(3) 

224  ER  a  TEMP  a  BY<2) 

RY(3)  «  Y  ♦  BY ( 1 ) 

IF  (SIGN)  225,226,226 

225  BY(3)  *  -BY(3) 

226  BY (3 )  ■ ( ( bR  -  BY(3)))*  Y 
IF  (FLX  -FLl)  222,228,228 

228  IF  (SIGN)  430,430,100 

222  BJ(D  a  8  J  (2  ) 

BY ( i )  «  BY (2 ) 

Bj(2 )  a  B J ( 3 ) 

8Y(2)  a  8Y (3 ) 

FLX  *  FLX  ♦  1, 

GB  TB  221 

300  LX  a  L  ’  1 
SEXIT  a  1, 

3X1  OB  320  I  a  1,  3 

301  SUM?  ■  1, 

SUM1  a  1. 

TfcRMl  *  1, 

TfcRM2  a  1, 

FM  a  1, 
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TL  a  2  *  LX 
3U2  TM  s  2,  *  FM 

TfcRMl  a-TbRHl  •  ZSG  /  ( TH  •  ( TU  ♦  TM  ♦  1.))  *  SEXIT 

TbF<M2  s  TERM2  *  7bQ  /  {TM  *  (TU  «  TM  ♦  1,))  *  SfcXl ▼ 

0LO1  *  SUMl 

0Ln2  s  SUM2 

SUMl  *  SUMl  ♦  TfcRMl 

SUM2  »  SUM2  ♦  TERM? 

DIFF  a  A8SFU0LD1  -  SUMl)  /  SUMl)  ♦  ABSFC(Bl02  -  SUM? ) 
IF  (OlFF  *  ,000001)  3u4,  306,  304 
•30  4  FM  sFM  ♦  1, 

G0  T0  302 

306  IXa  2  *  LX-  1 
ASM  s  -1, 

IF  < *-X )  310#  310,  307 

307  P@  30b  J  a  l,  IX,  ? 

AU  a  J 

308  ASM  s  ASM  «  AJ 

310  ZYL  *  (z  /  Y)  *•  LX 
ZYU=  ZYL  •  <Z  /  Y) 

IF  (LX)  312,  312,  314 
312  ASMl  a  1, 

G0  T0  316 
314  ALX  a  LX 

ASMl  a  -ASM  «  (2,  *  ALX  ♦  l,) 

316  0U<1)  a  ZYl  •  SUMl  /  ASMl 
BY(I)  a  ASM  •  SUM2  /  /YU 
320  LX  a  l,X  ♦  1 

IF  (SfcXIT)  430.430,900 
400  SIGN  a  -1, 

ER  a  fcXPF  <Z0) 

ERL  »  1,  /  ER 

SNH  a  ,5a  { ER  -  tRl.) 

CSW  a  ,5*  1ER  *  tRU 
B J ( 1 )  a  SNH  /  ZB 
RY ( 1 )  a  Y  •  CSH 
G0  TQ  210 
500  LX  a  L  -  l 
SEX  I T  a  -l, 

G0  T0  391 
430  KNT  a  2 

0Y ( 1)  a  Y/Z0 

BY ( 2 )  a  Y*Y* ( 1 , *ZB )  /  UB*Z0) 

5b3  BY ( 3 )  *  < 2 , *FL0ATF <KN T)  -1,)*Y*BY{2)  /  ZB  ♦  Y*Y*BYU> 
IF (KNT-L-1)  582,5ttl,5bl 
5tt?  BY  (1 )  a  BY  <  2 ) 

BY  C 2 >  a  RY<3> 

KNT  a  KNT  ♦  1 
G0  T0  583 

581  EXZ  a  EXPFI-ZB) 

00  584  N  <  1,3 
584  BY ( N )  a  BY ( N )  •  EXZ 

IF  (ZB-l.)  600,580,580 
5tt0  IF  <Z*4.aKUR>  700,600,600 
600  ER  a  2,*  FL  *  1, 

RJP  a  (FL  •  R J ( 1 )  /  Y  ♦  FLl  *  9J(3)  •  Y)  /  ER 
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QYP  aMFL  *  RYU)  •  Y  ♦  FLt  •  *Y<3)  /  Y)  /  E9 

Lac*i 

I ►  ( L  ,  EQ , 0  >  Gd  T 0  110U 

fljPP  a  2.  •  PJP  /  7  -  (1,  ♦  M.  *  FL1  /  ZSO)  •  dJ<?) 

RYP?  »  2.  *  HYP  /  7  -  Cl.  ♦  M.  «  f'Ll  /  Z40)  *  «Y(2) 

Gtf  T0  950 
POO  A  a  0,1 
IZE  a  U 
8  a  0,33 

LI  a  l,n 

FN1  a  L 
IDX  a  *1 

FN2  a  Z  -  0,5  ♦  SuRTF(3U,U  *  a  *  l) 

Zb  a  0, 

IJl  *  2,  •  7/ <2.*  Ml  ♦  1.) 

U2  a  2,  •  7/ (2 . •  F N2  ♦  1.) 

Ml  a  Ml  ♦  30,  •  (A  *  R*tji  •  ( ?  , -ut  *U1 )  /  (2  .  Ml.  -  U1  •  U1 ) )  ) 

M2  a  M2  ♦  30.  •  (A  ♦  R*U2  *  ( 2  , -U?*U2 )/ (2 .  Mi.  -  U2  •  U2 ) )  > 

AL  a  L 

IF  (F.N2-AU  801, 801, 8U? 

801  AM  a  Ml  ♦  l 
M  a  Ml 

Gfc)  T0  H03 

802  IF  (Ml  -  M2 }  801,8ul#«05 

805  AM  a  M2  *1 
M  a  M2 

803  R <  3 )  a  2/(2.  *  AM  *  3,  ) 

N  a  M 

800  AN  a  N 

R<2)  «  Z/(2,  *  AN  *  3,  -  Z  •  R(3)) 

N  a  N  -1 

IF  (L  -  N)  806,806,80/ 

806  R (3 )  a  R<2) 

Gfc)  T0  R08 

e07  BJ(3)  a  R<2) 

Bj<2>  a  1, 

AL  a  2*L  *1 
BJ(1)  a  AL/Z  -  H<2> 

809  LA  a  *L 
LAla-1,-1 

ALPHA  a  ZaZ*(aJ(2)*dY(l)*Y**LA-dJ(l)*dY(2)*Ya*LAl) 

NA  a  L-*l 
00  812  N3 1 • 3 

Bj(N)a(l./(Y**NA*ALPHA) )  *3J(N) 

512  NAaNAal 
G0  T9  900 

700  A  a  0,1 
8  a  ,25 

FN  a  « , 5  *  SQRTF  <30 ;  •  ti  *  ZB) 

IF  (FN-FL)  701.702.702 

701  FN  a  FL 

702  U  «  2,  •  ZB  / C 2 1  »'  FN  ♦  1.) 

M  a  FN  ♦  30,  ♦  (A  ♦  8  •  U)  ♦  1, 

R ( 3 )  a  ZB/  2.  •  FL0A  rF ( M )  ♦  1, 

M  a  M  -  1 
705  AN  a  M 
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R(2)  *  ZB/<2,  •  AN  ♦  1,  ♦  ZB  •  R<3)) 

M  *  M  -  1 

IF  (M  -  L  -1)  703(704,704 
704  R(3)  *  R ( 2 > 

G0  T0  705 
703  BJ<3)  *  R<2) 

BJ ( 2 )  >  1, 

8J(1)  • < 2 , •FU  ♦  1,)/ZB  *K<2) 

AUPHA«Z8*ZB*(BJ(2)*BY<1)/Y**U  *  9J<l)*BY(2)/Y**Ll ) 
NA*FL*1. 

38  712  Nsi , 3 

RU(M)«  (1,/(Y**NA*AUPHA)>  •  RJ(N) 

712  NA=NA*1 
G0  T0  600 

9UQ  ER  •  l,  /  <2,  *  FI  ♦  1.) 

BjP  a  ER  *  (FL  <  BJ  ( 1 )  /  Y  *  FU1  •  BJ(3)  •  Y) 

BYP  1  6R  ‘  (FI,  *  BY(l)  *  Y  »  F(.l  •  8Y{3)  /  Y) 

UBC*2 

IF  <U,EQ.Q)  GB  T0  1100 

BJPP  «  (FU  •  FL1  /  ZSU  -  1.)  •  BJ(2)  ,  2.  •  BjP  /  Z 

BYPP  «  (FU  •  KL1  /  ZSU  -  1.)  •  BY (2 )  -  2.  •  BYP  /  Z 

950  BA  a  BJ(2) 

BB  n  BY ( 2 ) 

1000  Y  a  YSAVE 
RETURN 

1100  IF  (IKQ.Lc, 1)1120(1110 
1110  IF  (QBO.UT, 0)1150, 1120 
1120  BAaSIN(Z)/Z 
RBs-CBS ( Z )  /Z 
IF  (UBU-DH30, 1130, 1140 
1130  BJPPa2.*BuP/Z-(l,*FU*FUl/ZSU)*BA 
BYPPa2,*8YP/Z-(l,*FU*FUl/ZSQ)*BB 
RETURN 

1140  8UPPa(FL*FLl/ZSQ-l, )«BA-2, *BJP/Z 
BYPPa(Fu*FUl/ZSQ«l, )*BB-2,*BYP/Z 
RETURN 

1150  9AaSlNH(Z)/Z 

Pla3, 1415926535 
9BaEXP(-Z)*PI/(2,*Z) 

IF  <183-1)1130,1130,1140 
END 


FUNCT 1 0N  SINH(Z) 

SINHa(EXP(Z)-fcXP(-7))/?, 

ENO 
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